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Abstract 

We introduce and discuss certain kinetic models of (continuous) opinion forma- 
tion involving both exchange of opinion between individual agents and diffusion of 
information. We show conditions which ensure that the kinetic model reaches non 
trivial stationary states in case of lack of diffusion in correspondence of some opinion 
CN| ' point. Analytical results are then obtained by considering a suitable asymptotic limit 

of the model yielding a Fokker-Planck equation for the distribution of opinion among 
' individuals. Numerical results on the kinetic model confirm the previous analysis. 
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1 Introduction 



Microscopic models of both social and political phenomena describing collective behaviors 
and self-organization in a society have been recently introduced and analyzed by several 
authors EEH EH HH HHj • The leading idea is that collective behaviors of a society 

composed by a sufficiently large number of individuals (agents) can be hopefully described 
using the laws of statistical mechanics as it happens in a physical system composed of 
many interacting particles. The details of the social interactions between agents then 
characterize the emerging statistical phenomena. 

Among others, the modelling of opinion formation attracted the interest of a increasing 
number of researchers (cfr. [Hj IT% 1 127j and the references therein). The starting point 
of a large part of these models, however, is represented by a cellular automata, where 
the lattice points are the agents, and where any of the agents of a community is initially 
associated with a random distribution of numbers, one of which is the opinion. Hence 
society is modelled as a graph, where each agent interacts with his neighborhoods in 
iterative way. 

Very recently, other attempts have been successfully applied p][22], with the aim to 
describe formation of opinion by means of mean fields model equations. These models 
are in general described by systems of ordinary differential equations or partial differential 
equations of diffusive type, that can in some case be treated analytically to give explicit 
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steady states. In 0, attention has been focused on two aspects of opinion formation, 
which in principle could be responsible of the formation of coherent structures. The first 
one is the remarkably simple compromise process, in which pairs of agents reach a fair 
compromise after exchanging opinions [21 El El EE H21 123 EH! • 

The second one is the diffusion process, which allows individual agents to change their 
opinions in a random diffusive fashion. While the compromise process has its basis on the 
human tendency to settle conflicts, diffusion accounts for the possibility that people may 
change opinion through a global access to information. In the present time, this aspect 
is gaining in importance due to the emerging of new possibilities (among them electronic 
mail and web navigation |2L)|). 

Following this line of thought, we consider here a class of kinetic models of opinion 
formation, based on two-body interactions involving both compromise and diffusion prop- 
erties in exchanges between individuals. Compromise and diffusion will be quantified by 
two parameters, which are mainly responsible of the behavior of the model, and allow for a 
rigorous asymptotic analysis. In consequence of our assumptions on the microscopic inter- 
action, in various relevant cases the model will satisfy mass and momentum conservation, 
which are the starting point for studying the asymptotic behavior. 

In this direction, we shall show that the kinetic model gives in a suitable asymp- 
totic limit (hereafter called quasi-invariant opinion limit) a partial differential equation 
of Fokker-Planck type for the distribution of opinion among individuals. Similar diffu- 
sion equations were obtained recently in [22] as the mean field limit of the Ochrombel 
simplification of the Sznajd model 27 . 

The equilibrium state of the Fokker-Planck equation can be computed explicitly and 
reveals formation of picks in correspondence to the points where diffusion is missing. 

The mathematical methods we use are close to those used in the context of kinetic 
theory of granular gases, where the limit procedure is known as quasi-elastic asymptotics 
|17| l2*%] . We mention here that a similar asymptotic analysis was performed on a kinetic 
model of a simple market economy with a constant growth mechanism I21j . showing 
formation of steady states with Pareto tails JH]- In this context, the mean field approx- 
imation leads to the same Fokker-Planck equation [S1E3], showing consistency between 
kinetic and stochastic approaches. 

The paper is organized as follows. In the next section we introduce the binary inter- 
action between agents, which is at the basis of the kinetic model. The main properties of 
the model are discussed in section 3. These properties justify the quasi-invariant opinion 
limit procedure, performed in section 4. The limit is illustrated by several choices of the 
diffusion in section 5. 

2 Description of the kinetic model 

The goal of the forthcoming kinetic model of opinion formation, is to describe the evolution 
of the distribution of opinions in a society by means of microscopic interactions among 
agents or individuals which exchange information. To fix ideas, we associate opinion with 
a variable w which varies continuously from —1 to 1, where —1 and 1 clearly denote the two 
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(extreme) opposite opinions. We will moreover assume that interactions do not destroy 
the bounds, which corresponds to impose that the extreme opinions can not be crossed. 
This crucial rule emphasizes the difference between the present social interactions, where 
not all outcomes are permitted, and the classical interactions between molecules, familiar 
to people working in kinetic theory of rarefied gases UJ. 

Let X = [— 1,+1] denote the interval of possible opinions. From a microscopic view 
point, we describe the binary interaction by the rules 

w 1 = w — r yP(\w\)(w — w*) + 7]D(\w\) 

(1) 

= -7P(|w*|)(w* -w) +77*D(|u;*|) 

where the pair (w,w*), with w,w* £l denotes the opinions of two arbitrary individuals 
before the interaction and (v/,w^) their opinions after exchanging information between 
them and with the exterior. In (^Q) we will not allow opinions to cross boundaries, and 
thus the interaction takes place only if both w', w'* G T. In JTJ) the coefficient 7 6 (0, 1/2) 
is a given constant, while 77 and 77* are random variables with the same distribution with 
variance a 2 and zero mean, taking values on a set B C 1R. The constant 7 and the vari- 
ance a 2 measure respectively the compromise propensity and the modification of opinion 
due to diffusion. Finally, the functions P(-) and D (•) describe the local relevance of the 
compromise and diffusion for a given opinion. 

Let us describe the details of the interaction in the right hand side. The first part is 
related to the compromise propensity of the agents, and the last contains the diffusion 
effects of external events. Note that the pre-interaction opinion w increases (getting closer 
to u>*) when > w and decreases in the opposite situation. The presence of both the 
functions P(-) and -D(-) is linked to the hypothesis that the availability to the change of 
opinion is linked to the opinion itself, and decreases as soon as one gets closer to extremal 
opinions. This corresponds to the realistic idea that extremal opinions are more difficult 
to change. We will present later on various realizations of these functions. In all cases, 
however, we assume that both P(|u;|) and |) are non increasing with respect to \w\, 
and in addition < P{\w\) < 1, < D(\w\) < 1. 

In absence of the diffusion contribution (77,77* = 0), Q implies 

w' + w'* = W + 7/J* + j(lV — W*) (P(\w\) — P(\w*\)) 

(2) 

v/-vi = (l-2 7 (P(\w\)+P(\w*\)))(w-w*). 

Thus, unless the function P(-) is assumed constant, P = 1, the total momentum is not 
conserved and it can increase or decrease depending on the opinions before the interaction. 
If P(-) is assumed constant, the conservation law is reminiscent of analogous conservations 
which take place in kinetic theory. In such a situation, thanks to the bounds on the 
coefficient 7, equations Q correspond to a granular gas like interaction (or to a traffic 
flow model ) where the stationary state is a Dirac delta centered in the average opinion 
(usually referred to as synchronized traffic state in traffic flow modelling). This behavior is 
a consequence of the fact that, in a single interaction, the compromise propensity implies 
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that the difference of opinion is diminishing, with \w' — w'J = (1 — 27) \w — w*\. Thus 
all agents will end up in the society with exactly the same opinion. Note that in this 
elementary case a constant part of the relative opinion is restituted after the interaction. 
This property does not remain true if the function P depends on the opinion variable. In 
this case 

\w' - wi\ = (1 - 2 7 (P(M) + P(K|))) \w-w*\. 
In fact, since 7 G (0, 1/2) and < P(\w\) < 1, < D(\w\) < 1, 

< e(w,w m ) = 1 - 2j(P(\w\) + P(K|)) < 1. 

Hence the general case corresponds to a granular gas interaction with a variable coefficient 
of restitution |28j . 

We remark moreover that, in absence of diffusion, the lateral bounds are not violated, 
since 

w' = (1 — jP(\w\))w + 7-P(|w|)u)* 

(3) 

wi = (l- 7 P(K|))w*+7P(K|)w 

implies 

max {\w'\, \w'J} < max{|u;|, . 

Let f(w, t) denote the distribution of opinion w G T at time £ > 0. A direct application 
of standard methods of kinetic theory of binary interactions 6 allows to recover the time 
evolution of / as a balance between bilinear gain and loss of opinion terms, described by 
the integro-differential equation of Boltzmann type 

% = J J X ^jf^fC™*) - 0/M/(wo) dw* drj d V *, (4) 

where ('w,' w*) are the pre-interaction opinions that generate the couple (w, w*) of opinions 
after the interaction. In Q J is the Jacobian of the transformation of (w, u;*) into (w' , w*), 
while the kernels ' (3 and (3 are related to the details of the binary interaction. 

As usual in classical kinetic theory of rarefied gases, the interaction integral on right- 
hand side of Q represents the instantaneous variation of the distribution of opinion, due 
to the binary exchanges of information. The presence of the Jacobian J, guarantees that 
equation © preserves the mass (total opinion), for any choice of the rate function (3. The 
transition rate is taken of the form 

/W.HKM) = ©fa)6fa*MM < !MKI < !)> ( 5 ) 

where x(.A) is the indicator function of the set A, and O(-) is a symmetric probability 
density with zero mean and variance a 2 . The rate function 0( w ,w r )-+(w',wi) characterizes 
the effects of external events on opinion through the distribution of the random variables 
and 0* and takes into account the hypothesis that bounds can not be violated. We 
remark that in principle the support B of the symmetric random variable is a subset of 
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X, to prevent diffusion to generate a complete change of opinion. This property can be 
weakened by assuming for example diffusion as a random variable normally distributed 
but well concentrated on zero. 

For a general probability density &(■), the rate function (5 depends on the opinion 
variables (w,w*) through the indicator functions \- This fact reminds a similar property 
of the classical Boltzmann equation jH) |§] , where the rate function depends on the relative 
velocity. As we shall see, a simplified situation occurs when a suitable choice of the function 
D{-) in Q coupled with a small support B of random variables implies that both \w'\ < 1 
and \w+\ < 1, and the kernel (3 does not depend on the opinion variables (w,w*). In this 
case the kinetic equation © is the corresponding of the classical Boltzmann equation for 
Maxwell molecules jl], which presents several mathematical simplifications. In all cases 
however, methods borrowed from kinetic theory of rarefied gas can be used to study the 
evolution of the function /. 



3 Simplifications and main properties of the model 

The main problem in opinion dynamics is the formation of stationary profiles for the 
opinion. In the kinetic picture this corresponds to an investigation of the large time 
behavior of the density of opinion f(w, t). To investigate in detail the large-time behavior, 
a preliminary analysis of equation (jlj) is needed. We will start this analysis by introducing 
some notations and by discussing the main properties of the kinetic equation. 
Let Q(f, f) denote the interaction integral, 

Q(fJ)(w) = ^J^hjf('w)f('w,)-f3f(w)f(w^dw,di 1 dr,,. (6) 

Let A4q(A) the space of all probability measures taking values in A C M and by 



M P (A) = |© e Mo : J \w\ p dQ(w) < +oo, P > j , 



(7) 



the space of all Borel probability measures of finite momentum of order p, equipped with 
the topology of the weak convergence of the measures. 

Let Tsil), be the class of all real functions feonl such that h(±l) = h'(±l) = 0, and 
h,( m \v) is Holder continuous of order 5, 

\\h [m) \\ s = sup- Y \/j < °°> 8 ) 

V -L W \V-W\ a 

the integer m and the number < 8 < 1 are such that m + 5 = s, and denotes the 
m-th derivative of h. 

In the rest of the paper we will assume that the symmetric probability density 0(7/) 
which characterizes the diffusion of information belongs to A^+aj for some a > 0. More- 
over, to simplify computations, we assume that this density is obtained from a given 
random variable Y with zero mean and unit variance, that belongs to M.2+ a - Thus, O of 
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variance a 2 is the density of oY . By this assumption, we can easily obtain the dependence 
on a of the moments of O. In fact, for any p > such that the p-th moment of Y exists, 



f \r 1 \P@(r ] )dr ] = E{\aY\ p )=a p E{\Y\ p ). 



By a weak solution of the initial value problem for equation (flf). corresponding to the 
initial probability density fo(w) G Mq(I), we shall mean any probability density / G 
C 1 (M+, .Mo (2")) satisfying the weak form of the equation 



dl 



cj ) (w)f(w,t)dw = (Q(f,f),0) = 
(3( w ,w*)^(w>, w i)f(w)f(w*)((/)(w') - 4>(w))dw«dwdr]dri*, (9) 



1 1 2 Jb 2 

for i > and all 4> G J- p (2), and such that for all cj) G Fp(I) 



lim / (f>(w)f(w,t)dw= / (p(w) fo(w) dw . (10) 

The form Q is easier to handle, and it is the starting point to explore the evolution of 
macroscopic quantities (moments). By symmetry reasons, we can alternatively use the 
symmetric form 

(11) 

(4>(w) + 4>{ w '*) — 4>( w ) — 0(^*))^w* dwd-qdrj*. 

Existence of a weak solution to the initial value problem for equation (j3J) can be easily 
obtained by using methods first applied to the Boltzmann equation 6 . On the other hand, 
for a general kernel it appears extremely difficult to describe in detail the large-time 
behavior of the solution. 

For this reason, at first we restrict our analysis to the cases in which the kernel (5 
does not depend on the opinion variables (Maxwellian case). In this direction, let us 
briefly discuss the importance of the support B of the probability density function 0(7/) 
in connection with the possible simplification of the kernel. To clarify the point, let us set 

D(\w\) = 1 - \w\. 

This function satisfies all the requirements we fixed in the previous Section. Then, since 

(i - yp(\w\))w + 7P(M)™* + (i - H)v < (i - 7-P(M))«> + iP(H) + (i - H)v, 

in order that \w'\ < 1 it suffices that 

(1 - jP(\w\))w + jP{\w\) + (1 - \w\)t] < 1, (12) 
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or, what is the same 

(i-M>/<(i-7^(M))(i-«0. (13) 

Since P(|u;|) < 1, bound is verified for all w > 0, as soon as r) < 1 — 7. Analogous 
result holds if ui < 0. Hence, if D(|u;|) = 1 — \w\ and B = (—(1 — 7), 1 — 7), both w' and 
w'x belong to the right interval. 

Remark 3.1 Any choice oi D(\w\) and B which are suitable to preserve the lateral bounds 
of extreme opinions allows to study in detail the dynamics of the model with a significant 
simplification. In these cases, the kernel (3 defined in © simplifies to 

/3(te,w*)-K,0 = PiViV*) = ®(v)®{v*)- 

In the rest of the paper we will limit ourselves to such type of kernels. As briefly explained 
in the previous Section, this assumption is the analogue to Maxwell molecules interaction 
in the Boltzmann equation jSJ. 

From © (or equivalently from Q11JI ) conservation of the total opinion is obtained for 
(p(w) = 1, which represents in general the only conservation property satisfied by the 
system. The choice <fi(w) = w is of particular interest since it gives the time evolution of 
the average opinion. We have 

— / wf(w,t)dw= I [ /3(ri,i7*)f(w)f(w*)^(P(\w\)w* - P(\w\)w)dw*dwdridr]* 
dt Jx Jiz Jb 2 



+ / P(v,V*)f( w )f( w *)' l l D (\ w \)dw*dwdr]dri* 
Jx 2 Jb 2 

The first integral on the right-hand side represents the contribution of the exchange of in- 
formation to the variation of momentum. In case P(|u>|) = 1, this contribution disappears, 
since, by symmetry, 

/ V*) f ( w ) f ( w *)li w * — w)dw*dwdr] dr]* = 0. 

x 2 Jb 2 



In this case 

d_ 

di 



wf(w, t) dw 



x 



= / / v®(v)®(ri*M\w'\ < 1MKI < l)D(\w\)f(w)f(w,)dw*dwdn dr/* = 0. (14) 
Jx 2 Jb 2 

since the mean value of the density B is zero. This shows that P constant implies that 
the average opinion is conserved. The situation changes when P is not constant. In this 
case, the time evolution of the average opinion is given by 



d 
dt 



w f (w , t) dw = 7 / P(\w\)f(w) dw / wf(w)dw — r y / wP(\w\)f(w) dw. (15) 
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Note that equation (fT5|) is not closed. 
Let us fix now cj>(w) = w 2 . We have 

— f w 2 f(w,t)dw = - f f Q(r])Q(7]*)f(w)f(w*)(w' z + w / ^ — w 2 — wl)dw*dwdridri* 

dt Jj 2 Jj2 J B 2 

(16) 

Taking in mind that O has zero mean and variance u 2 , by easy computations one shows 
that 

- / / @(rj)Q(r]*)f(w)f(w*)(w' 2 +w'J 2 -w 2 -wl)dw*dwdr]dri* = 



i 2 jb 2 

2 I r>n„..\\2f„.. „,. \2 . 



! P{\ 

Jj 2 



7 / P(\w\) (w — w*) f(w)f(w*)dwdw 1r 



2 7 J ^ P(\w\)w(w - w*)f(w)f(w*) dwdw, + a 2 J D(\w\) 2 f(w)dw. (17) 
The choice -P(|u>|) = 1 leads to the simpler evolution equation 



d_ 

dt 



w 2 f(w) dw — m 2 



w f(w, t) dw = —27(1 — 7) 
where m is the constant value of the average opinion 

m = wf(w, t)dw. 



+ a 2 j^D(\w\) 2 f(w)dw, (18) 



J% 

Since \w\ < 1, the boundedness of the mass implies that all moments are bounded. This 
implies that, in all cases, we can draw conclusions on the large-time convergence of the 
class of probability densities {f(w,t)} t>0 . By virtue of Prokhorov theorem (cfr. |TI]) the 
existence of a uniform bound on moments implies in fact that this class is tight, so that 
any sequence {f(w, t n )} n>0 contains an infinite subsequence which converges weakly to 
some probability measure f^. 



4 The quasi-invariant opinion limit 

The analysis of the previous Section shows that in general it is quite difficult both to study 
in detail the evolution of the opinion density, and to describe its asymptotic behavior. For 
a general kernel one has in addition to take into account that the mean opinion is varying in 
time. As is usual in kinetic theory, however, particular asymptotics of the equation result 
in simplified models (generally of Fokker-Planck type), for which it is relatively easier to 
find steady states, and to prove their stability. These asymptotics are particularly relevant 
in case they are able to describe with a good approximation the stationary profiles of the 
kinetic equation. In order to give a physical basis to these asymptotics, let us discuss 
the interaction rule Q from a slightly different point of view. For the moment, we will 
assume P(|w|) = 1, so that conservation both of mass and momentum holds. The case of 
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a general will be treated subsequently. Let us denote by E(X) the mathematical 

expectation of the random variable X. Then the following properties follow from |T)l 

E[w' + w' t ] = w + w„ %'-t»;] = (l-27)(w-4 (19) 

The first equality in (|19|) describes the property of mean conservation of opinion. The 
second refers to the compromise propensity, which plays in favor of the decrease (in mean) 
of the distance of opinions after the interaction. This tendency is a universal consequence 
of the rule (^Q), in that it holds whatever distribution one assigns to 0, namely to the 
random variable which accounts for the effects of the external word in opinion formation. 

The second property in IJ19JI is analogous to the similar one that holds in a collision 
between molecules in a granular gas. There the quantity e = 27 is called "coefficient of 
restitution", and describes the peculiar fact that energy is dissipated |28] . 

We consider now the situation in which most of the interactions produce a very small 
exchange of opinion (7 — ► 0), while at the same time both properties (|19ft remain at a 
macroscopic level. This corresponds to pretend that, while 7 — > 0, 

w + w*) f (w) f (w*)dwdw* = 2 / wf(w)dw = 2m(t) (20) 
x 2 ' Ji 

remains constant, and 

1 

2 



w — w*) f(w)f{w*)dwdw* = / w f (w) dw — m = C f (t) (21) 
x 2 Ji 



-2 7 (1 - 7 )C > (t) + a 2 J r D(\w\ff(w)dw. (22) 



varies with time, and decays to zero when the diffusion is not present (i.e. a = 0). 

Since in our case the kernel (3 does not depend on the opinion variables, (j!4j) implies 
that m(t) = rriQ independently of the value of 7. Moreover, using the computations of the 
previous Section, one obtains that Cf(t) varies with law 

dC f (t) 
dt 

Hence, if we set 

r = 7 t, 9 (w,t) =f(w, t), (23) 
which implies fo(w) = go(w), it follows 

d £M =-2(1- 7) C g {r) + - f D(\w\) 2 f(w)dw. (24) 
dr 7 Ji 

Letting now both 7 — > and a — > in such a way that cr 2 /7 = A, (|24[) becomes in the 
limit 

dC gi.V _ -2C 9 {t) + A f D(\w\) 2 f(w)dw. (25) 



1 



dr gy 

This argument shows that the value of the ratio cr 2 /7 is of paramount importance to 
get asymptotics which maintain memory of the microscopic interactions. It is remarkable 
that, thanks to (|2.3|) . t = r/7, so that the limit 7^0 describes the large-time behavior 
of f(v,t). On the other hand, since f(w,t) = g(w,r) the large-time behavior of f(w,t) is 
close to the large-time behavior of g(w,r). 
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Remark 4.1 The balance 7 — ► and a — > in such a way that <7 2 /7 = A, allows to 
recover in the limit the contributions due both to compromise propensity and the diffusion. 
Other limits can be considered, which are diffusion dominated (c 2 /7 = 00) or compromise 
dominated (c 2 /7 = 0). As we shall present in the next section, however, the formation of 
an asymptotic profile for the opinion is linked to the first balance PP. 

In the remainder of this section, we shall present the rigorous derivation of a Fokker- 
Planck model, starting from the Boltzmann equation for the opinion density g(w, r), when 
both 7^0 and a —* in such a way that <r 2 /7 — ► A. For the sake of simplicity, we will 
assume that f (|uu|) = 1. This type of analysis is close to the one described in |7j for a 
kinetic model of wealth distribution in a open economy. 

The scaled density g(v,r) = f(v,t) satisfies the equation (in weak form) 

-7- [ g{w)cj)(w) dw = - f f @(rj)@(ri*)g(w)g(w*)(<j)(w r ) — (f>(w))dw*dwdr] drj*. (26) 

dT J J 7 Jj2 J B 2 

Given < 5 < a, let us set (j) E (I) ■ 
By©, 

w — w = 7(w* — w) + 7jD(\w\). 
Then, if we use a second order Taylor expansion of (f> around w 

<j)(w') - (j)(w) = (7(10, -w) + i]D(\w\))<t> f (w) + - (7(10* -w) + r]D{\w\)) 2 <j>"(w), 
where, for some < 6 < 1 

w = Qw + (1 — 6)w. 
Inserting this expansion in the collision operator, we get 

4- I g(w)(j)(w)dw = - I I e(r})@{r]*){( 1 (w*-w)+r}D(\w\))<t)\w)+ 
dr J j 7 J T 2 J B 2 

+^(l( w * ~ w ) + r]D(\w\)) 2 (p'\w)]g{w)g(w*)dw*dw d-qdri* + R('y,<r), (27) 

where 

R{rf,G) = ± f f 9( r] )e(r ] *)( 1 (w*-w)+ V D(\w\)) 2 - 
27 Jt? Jb* 

■ (4>"{w) — 4>"{w)) g(w)g(w*)dw* dw dr] drj* 
Since (j) E J r 2+s{Z), an d \w — w\ = 6\w' — w\ 

\4>"{w) - <P"{w)\ < U"\\s\w - w\ s < U"\\s\w' - w\ s . (28) 

Hence 

\R(%a)\ < M]k f f 9(^9^) • 
27 Jt 2 Jb 2 

•|7(tu* — w) + r]D(\w\)\ 2+s g(w)g(w*)dw* dwdr/dr]* 
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By virtue of the inequality 



| 7 K -w)+ r)D(\w\)\ 2+s < 2 1+s (| 7 (w* " ^)| 2+5 + \vD(\w\)\ 
<2 2+s 1 1+s + 2 1+S \r ] \ 2+S , 

we finally obtain the bound 



2+5 



< 



\R( 7 ,a)\<2 1+5 U"\\s{l 1+5 + 



1 

27 



rj\ 2+5 0(rj)dn 



(29) 



Since is a probability density with zero mean and A7 variance, and @ belongs to A^2+«, 
for a > 5, 



L 



\?]\ 2+5 e(r])d7] = E 



A 7 y 



2+8 



(A 7 ) 1+5/2 £ (\Y\ 2+S ) 



and E \ \Y^ + j is bounded. Using this equality into (|29|) one shows that i?( 7 , a) converges 

to zero as as both 7 and a converge to zero, in such a way that a 2 = A 7 . Within the same 
scaling, 

e(r))@(r]*)[(j(w*-w)+riD(\w\))<l>'(w)+ 



lim — 

7^0 7 Jj2 J B 2 



1 



+-( 7 (w* — to) + r]D(\w\)) <p"(w)]g{w)g{w^)dw tf dwdr]dri^ 



(m - + -D(\w\) 2 (f>"(w) 



g{w)dw 



(30) 



Considering that cp £ J- S (T), we can integrate back by parts. This shows that the right- 
hand side of (|3U|) coincides with the weak form of the Fokker-Planck equation 



dr 2 dw 2 



'D{\w\) 2 g) + ^-{{w-m)g). 



(31) 



Last, since the solution to the kinetic model conserves mass and momentum, while the 
second moment is uniformly bounded in time, conservation of both mass and momentum 
pass to the limit. We remark that these conservations are difficult to prove directly on the 
Fokker-Planck equation, due to the fact that 4>{v) = v does not belong to T S ^L). Hence 
we proved 



Theorem 4.2 Let the probability density fo £ and let the symmetric random 

variable Y which characterizes the kernel have a density in M.2+a, with a > 5. Then, 
as 7 — ► 0, a — > in such a way that a 2 = A 7 the weak solution to the Boltzmann 
equation for the scaled density g 7 (v,r) = f(v,t), with r = 7 i converges, up to extraction 
of a subsequence, to a probability density g(w,r). This density is a weak solution of the 
Fokker-Planck equation \31\) . and it is such that the average opinion is conserved. 
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5 Other Fokker-Planck models of opinion formation 

Theorem 14. 21 can be generalized in many ways. Always remaining with the simplification 
of Remark 13.11 we can choose a general function P(|u;|) into the interaction rule (fT]). The 
main difference with respect to the proof of Theorem l4.2l is the evaluation of the first order 
term into (|27|). which, using the fact that the mean value of is zero, reads 

P(\w\)(w* — w)<p' (w)g{w)g{w*)dw if dw = / P(\w\)(m(T) — w)cj)' (w)g(w) dw, (32) 

J 2 Jl 

where m(r) is the value of the average opinion at time r > 0, 

m(r) = J wg(w,r)dw. 

Note that, by (J23J) 

m(r) = J wg(w,T)dw = J wf(w,t)dw. (33) 

Hence by (|15|) the evolution of m(r) obeys the law 
dm(r) 



dr 



l Jl 



m(r) / P(\w\)g(w, r) dw — / wP(\w\)g(w, r) dw. (34) 



Finally, in the limit 7 — > we obtain that g(w, r) satisfies the Fokker-Planck equation 

% = \^ W w \) 2 9) + ^ (P(H)(™ ~ m(t))g) . (35) 



Remark 5.1 The presence of a general propensity function P(|u>|) introduces a difficult 
to treat nonlinearity into the Fokker-Planck equation. The nonlinearity is due to the fact 
that the average opinion is no more constant, and the evaluation of the drift term requires 
the evaluation of (|3*l)l . 

While to our knowledge the Fokker-Planck equation ()35|) has never been considered 
before, linked pure diffusion and drift equations have been recently introduced in |22j . 
These equations, in our picture, refer to diffusion dominated (<r 2 /7 = 00) or compromise 
dominated (o" 2 /7 = 0) limits. Looking at the proof of Theorem ()4.2|) . it is almost immedi- 
ate to conclude that the diffusion dominated limit takes into account only the second-order 
term into the Taylor expansion. To verify this, suppose that 




a < 1. 



Then we can set 

T = 7 Q i, 9 (w,t) = f(w,t), (36) 
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where now g(w, r) satisfies 

— / g{w)(f>{w)dw = \ [ ( e(r])e(r)*)[(j(w*-w)+vD(\w\))<t>'(w)+ 
dr Jx 7 Jx 2 JB 2 

+ 7^(l(. w * ~ w ) +r)D(\w\)) 2 (p'\w)]g{w)g(w*)dw Sf dwdr)dri* + R(j,a), (37) 

with obvious meaning of the remainder. Since a < 1, the first order term in the Taylor 
expansion vanishes in the limit, and g satisfies the diffusion equation 

| = W(«). (38) 

The choice 



D(\w\) = Vl-w 2 , A = 2, (39) 
brings to the diffusion equation 

dg d 2 



dr dw 2 



[(1 - w 2 )g] . (40) 



This diffusion equation has been derived in a mean field approximation 22 to describe 
the evolution of the Sznajd model in Ochrombel simplification |18j on a complete graph 
of iV nodes in case of two opinions. The same equation has been shown to describe the 
former model in case of a large number q of opinions. The variable now represents the 
mean distribution of occupation numbers and the limit is taken as both N and q tend to 
infinity at the same rate. 

Likewise, the compromise dominated (cr 2 /7 = 0) limit can be considered. In this case, 

a 2 



A, a > 1, 
and we can set 



r = 7i, g(w,r) = f(w,t). (41) 
The diffusion part disappears in the limit and we obtain the pure drift equation 

^ = JL(P(\ w \)(w-m(t))g). (42) 

Note that, due to (|41|) the evolution of the mean opinion m{t) obeys the law (|33|) . The 
choice 

P(\w\) = l-w 2 
has been considered in [21]. In this case 

f r = ^((l-w 2 )( W -m(t))g), (43) 

where 

— P^— = — m(r) / w 2 g(w,T)dw+ [ w 3 g(w,T)dw. (44) 
dr Jx Jx 
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Thus, our equation differs from the pure drift in magnetization obtained in |22| as the 
mean field limit of the Sznajd model [27] in case of two opinions. There the first-order 
partial differential equation reads 

| — £ ((1 (45) 

Note that the sign in front of the drift (|45|) is now opposite to the sign in (|43|) . and the 
equation is now linear in g, even if the evolution of the mean opinion is not closed. 



6 Stationary solutions of the Fokker-Planck opinion model 

In this Section we analyze in some details various cases of the interaction dynamics in 
the Boltzmann equation from which it is possible to derive a Fokker-Planck equation with 
a explicitly computable steady state. The structure of the steady state then represents 
the formation of opinion consequent to the choice of the interaction dynamics. In most 
cases, we are forced to suppose P(|u;|) = 1, which implies conservation of the average 
opinion, and Fokker-Planck (jUlj) as underlying quasi-invariant opinion limit. For any of 
these choices, we briefly discuss the link between D(|u>|) and the maximal support B of 
the diffusion variable. 

The first model of the diffusion dependence on opinion we propose is 

D(\w\) = l-w 2 . 

Since P(|k;|) = 1, the interaction rules are 

w' = w — 7(u> — w*) + 7/(1 — w 2 ) 

w'^ = w*- 7(w* - w) + ?7*(1 - w 2 ) 

In this case + \w\) < 1 — 7 implies \w'\ < 1. Hence, for a given opinion B = 
(— (1 — 7)/(l + \w\)i (1 — 7)/(l + l^l))) which shows that when the opinion w is close to 
zero and 7 is small, the effects of the global access to information can move the opinion 
towards extremals. This possibility reduces as soon as increases. The steady state 
distribution of opinion is a solution to 

±JL((l- w 2fg) + (w-m)g = (46) 

where m is a given constant (the average initial opinion) — 1 < m < 1. The solution to 
(|4l?j) is easily found as 

gooW = c m , x (1 + w )- 2+m/(2A) (1 + un,r W(2A) exp {"^=y } , (47) 

where the constant c mt \ is such that the mass of g^ is equal to one. Note that the presence 
of the exponential assures that ^(±1) = 0. The solution is regular, but not symmetric 
unless m = 0. Hence, the initial opinion distribution reflects on the steady state through 
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the mean opinion. In any case, the stationary distribution has two picks (on the right and 
on the left of zero) with intensity depending on A. 
A similar result is expected from the choice 

D(\w\) = 1 - \w\. 

As discussed in Section |2] \r/\ < 1 — 7 implies \w'\ < 1, and B = (—(1 — 7), 1 — 7). Once 
more, the support of the random variable, for 7 small, covers the whole domain of opinions. 
The steady state distribution of opinion is a solution to 

±JL((l-\w\) 2 g) + (w-m)g = (48) 

where m is a given constant (the average initial opinion) — 1 < m < 1. The solution to 
(|4*oT) is easily found as 



1 \ n 1 n~2-2/A / 1 - mw/\w\ 

goo(w) = c m , A (1 - M) p \ ~ 2A(1- H) 



(49) 



where, as usual, the constant c m) \ is such that the mass of g^ is equal to one. We remark 
that the low regularity of D(|w|) reflects on the steady solution, which has a jump in 
w = 0. The jump disappears only when the mean m = 0, and only in this case we have a 
symmetric distribution. As in the first case, the presence of the exponential assures that 
goo(±l) = 0, and the initial opinion distribution reflects on the steady state through the 
mean opinion. 

Last, we consider 

D(\w\) = \/l-w 2 

in the Fokker-Planck model The steady state distribution of opinion solves 

~ ((1 " w 2 )g) + (w - m)g = 0, (50) 

and equals 

^ \ l-(l+m)/A / 1 \ l-(l-m)/A 



ffooM = C m , A — • (51) 

\l + w J \1 — w J 

As usual, the constant c m) \ is such that the mass of is equal to one. Since — 1 < m < 1, 
#00 is integrable on I. Differently from the previous cases, however goo(w) tends to infinity 
as w — > ±1, and it has no peaks inside the interval I. The explanation comes out from a 
deep insight into the connection between D(|i/7|) and the support B in this case. 

It is immediate to verify that, in order to satisfy the constraint in (^Q), one can not 
choose directly D(|u)|) = Vl — w 2 . Within this choice, in fact, choosing w* = 1, the first 
equality in (^Q) gives for rj the upper bound 

.<d-7) {1 ~ W) 



y/l-W 2 ' 
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Since the right-hand side converges to zero as w — > 1, it follows that there is no way 
to satisfy the constraint. A different choice which gives in the limit the Fokker-Planck 
equation with = vl — w 2 , is the following. For any given 7, we set 

D(\w\) = ^(l-(l + 7 p K)+, P>0, 

where /+ denotes as usual the positive part of /, that is /+ = / if / > 0, while /+ = if 
/ < 0. In this case, one can show that it is sufficient for rj to satisfy the condition 

H< ai = ^=^'\ (52) 

to respect the constraint on post-interaction opinions. To give an example, let us assume 
that G is uniformly distributed on the interval — a 7 ,a 7 . Then a 2 behaves like ^f)/ 2 , and 
it is enough to set p = 2/3 to obtain A = 1. The previous discussion shows that the choice 
-D(|w;|) = \f\ — w 2 in the Fokker-Planck equation (|31|) corresponds to a kinetic interaction 
in which diffusion is of the order of j p , where p is taken so that <7 2 /7 tends to a finite limit 
A as 7 — > 0. In this case, the smallness of the interval of diffusion produces the peaks on 
w = ±1. 

A interesting feature of the Fokker-Planck equation 

| = w(p-^) + k«"-"W- (53) 

is that it leads to close evolution of moments. We remark that equation Q53|) is to be 
studied with the conservation of both mass and momentum. 



7 Conclusions 

We introduced and discussed here some kinetic models of opinion formation based on bi- 
nary interactions involving both compromise and diffusion properties in exchanges between 
individuals. A suitable scaling of compromise and diffusion allows to derive Fokker-Planck 
equations for which it is easy to recover the stationary distribution of opinion. Among 
these Fokker-Planck equations, one is emerging (|53[) and takes the role of the analogous 
one obtained in 0E] for the evolution of wealth. The main feature of this equation is that 
moments can be evaluated in closed form. Further numerical studies are in progress to 
understand the evolution of opinion density for various choices of the underlying functions 
P(|u>|) and In particular, when P(|u;|) is not linear, the evolution of moments is 

far from being completely understood. 
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